The independent scattering approximation in the saturated regime 
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We show that a saturable single-frequency elastic T-matrix approach to scattering of light by 
atoms agrees remarkably well with a master equation description in the regime of unsaturated 
atoms, or for large separation between the atoms. If the atoms are in each other's near-field and 
the saturation of the atoms is high enough, the two approaches yield different results. 
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I. INTRODUCTION 



Applications of multiple scattering of light by point- 
scatterers are often based on the assumption that the 
scatterers under consideration scatter light indepen- 
dently from one another. All the scattering properties of 
each individual scatterer can then be expressed by a sin- 
gle mathematical object called the scatterer's T-matrix 
0, 0] • This single-scatterer T-matrix is used as a build- 
ing block in multiple-scattering theories which can 
accurately describe many interesting phenomena 0, Q . 
However, the use of T-matrices usually implies that one 
considers only elastic scattering. Although the restriction 
of scattering events to the pure elastic ones is justifiable 
in some cases, it is not necessarily correct for, e.g., scat- 
tering of high-intensity incident light fields . 

In this paper, we will consider scattering of light by 
two atoms. The "atoms" could be implemented as any 
type of sub-wavelength quantum objects, for example: 
trapped atoms, quantum dots Q, trapped ions Q or 
dye molecules. If we model the atoms as point-dipoles, 
we can use the T-matrix formalism to describe scatter- 
ing of incident light. Alternatively, we can model the 
atoms as two-level systems, in which case multiple scat- 
tering of light can be studied using a master equation 
LLl| • The system is then characterized by a density 
matrix, representing the coherences and populations of 
the atomic levels. Scattering of incident light then deter- 
mines the evolution of the density matrix, while taking 
into account all orders of multiple scattering as well as 
inelastic scattering of light. 

The goal of this paper is to compare the scattering 
properties of two atoms as described by a simple T- 
matrix approach and a (computationally much more in- 
volved) master equation approach. We will show that 
a saturable single-frequency elastic T-matrix approach 
agrees remarkably well with a master equation descrip- 
tion in the regime of unsaturated atoms or for large sep- 
aration between the atoms. 
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FIG. 1: The system under consideration: two identical atoms 
1 and 2 are positioned in free space at a distance m from 
each other. Both atoms have a two-level internal structure. 
The upper level has a lifetime T^ a . Both atoms interact with 
an incident field; this interaction is depicted by the dashed 
arrows. 



The new aspect of this work is the comparison of two 
approaches which seem completely different at first but 
are in fact closely connected. We stress that the two ap- 
proaches themselves, reviewed for clarity in sections II 
and III of our paper jhave been extensively studied pre- 
viously (see, e.g., 

This paper is organized as follows. We start in Sec- 
tion II by describing the system we are considering. A 
master equation approach allows us to calculate how a 
near-resonant monochromatic incident field is scattered. 
In Section III, we consider the same system, and use a T- 
matrix approach to determine the scattering properties 
of the system. In Section IV we compare both methods. 
In Section V, finally, we discuss the applicability of our 
results to other atomic systems. 



II. THE DENSITY MATRIX APPROACH TO 
SCATTERING OF LIGHT 

We start by considering the system shown in Figure 
1. Two identical atoms are positioned in free space, 
separated by a distance r\i — \r\i\. Both atoms are 
modelled as a two-level system with upper level b and 
lower level a, separated by an energy difference hwba- 
Decay between both levels occurs according to the decay 



rate Tb a . A transition dipole moment d ab , i £ {1,2} is 
associated with both b — > a transitions. Both transition 
dipole moments are chosen to be equal in magnitude 
and orientation: 

d? ab = d a bn, Ia»| = 1, *e{l>2}. (1) 

An incident monochromatic field Eq with wave vector 
fco interacts with both atoms. The frequency loq = \ko\c 
of the incident field is tunable and chosen near the 
b — > a resonance (c is the velocity of light in free space) . 
For simplicity, we have taken the wave vector fco of 
the incident field perpendicular to r± 2 . This choice 
of fco ensures that both atoms always feel the same 
phase of the incident field, which makes the calculations 
somewhat easier without affecting the generality of our 
results. 

We will now derive the master equation of the system; 
the following derivation closely follows The total 

Hamiltonian Hq describing the energies of the systems, 
the electromagnetic field and interactions is, in the stan- 
dard electric dipole and rotating wave approximation, 
given by 

Ho = Hd + Hp + Hul + He>f- (2) 

The Hamiltonian of both individual atoms is given by 

H D = tiw ba (S+Si + S+S 2 ), (3) 

with Sf 1 the dipole raising and lowering operators of 
atom i. Both atoms are coupled to the three-dimensional 
multimode electromagnetic field with Hamiltonian 

Hp = {a{ x a kx + , (4) 

fcA 

where the operators aj^ and d^x respectively create and 
annihilate a photon in de mode (fc, A). The modes of the 
electromagnetic field are taken to be in a vacuum state. 
Both atoms also interact through the Hamiltonian 

H DL = hn((S+ + S+)e iUot + H.c.) , (5) 

with an incident field Eq with frequency luq, where the 
Rabi frequency is given by £1 = \d a \, ■ E[)\/h. Finally, 
the atoms interact with the multimode vacuum field 
through the Hamiltonian 

H DF = J2(» ■ (9kx(r!)S+ + g kx (r 2 )S+) +n.c}j, (6) 



where the mode function g k \ is given by 

with L 3 the quantization volume, Eq the vacuum permit- 
tivity and £k\ the unit polarization vector of the field 
mode (fc, A). 

The dynamics of the total system can be expressed 
in terms of the Master equation for the density matrix 
a. The Master equation is, in the standard Born and 
Markov approximation, written in the Lindblad form, 
given by [HEl 



— C n d<J + CdCf. (8) 

The non-dissipative part of the Lindblad operator can 
be written as 

C nd a = ~[H c +H ,a], (9) 
n 

with coupling Hamiltonian 

H c = h5^ (S+Sz + S+S{) , (10) 

while the dissipative part of the Master equation is given 
by 

C-do = — (jS^Si o — S\ <t S>i + S 2 S 2 <y — S 2 &S 2 ^J 

r( 12 ) / - - - - ~ \ 
— (sfS 2 & + &SlS 2 - 2S 2 &S+) + H.c. 

(11) 

It is clear from equations (|10|l and (|ll|l that the vacuum 
induces coupling between both atoms 0, 0, 0, Il7j . 
The parameters <5( 12 ' and r^ 12 ) are not associated with 
individual systems, but with the total system as a whole. 
The cross-damping r^ 12 ) represents incoherent coupling 
between both atoms through spontaneous emission 

r(12) = % E( d ' 9k\(ri)) (d ■ g h x(r 2 ))s(uj ba - «). 

fcA 

(12) 

The parameter <5( 12 ) represents coherent coupling 



through the vacuum, expressed by the frequency shift 



5< 



12) 



4-7T 



r+oo , ,3 

/ V (12) ( 



LO - LOba 



-)duj, 

(13) 



where V stands for Cauchy's Principal Value of the inte- 
gral. Using standard evaluating techniques and replacing 



E -> E ^3 / fc2dfc / d0 y Mew ( 14 ) 



we find, in accordance [2^ with 0] 

^(12) _ f__ CQS f x \ _|_ J_ gjnj'jj.) + _L cosfx)^ , 
4 V a; or a; 01 / 

T (12) = ^r ba (- wx{x) + ^ cos(x) - -1 sin(x)) , (15) 
2 Vie x or / 



where we defined the dimcnsionlcss parameter 
x = LL>ba,ri2/c, and where we used the fact that 
/Lt_Lri2. Both couplings and I^ 12 ) depend strongly 
on the distance between the atoms and the orientation 
of fi and J"i2. In Figure |3 we show how <5^ 12 ' and r^ 12 ' 
evolve as a function of x. For large separations (i» 1), 
the two couplings S^ 12 ^ and T^ 12 ) vanish, independent 
of the orientation of the dipole moments. However, 
for small separations [x <C 1), T^ 12 ) and S^ 12 ^ exhibit 
a different behavior: while S 1 -- 12 ^ diverges for small 
separations, T^ 12 ) reduces to the single-atom decay rate 

r. 

The standard method of solving the Master equation 
© is to calculate equations of motion for the density 
matrix elements and solve them by direct integration 
[H 03 . This procedure means a set of 16 linear 
differential equations has to be solved. Using the 
resulting density matrix, we will calculate the total 
field, resulting from the scattering of Eq. We therefore 
need to establish a relation between the total field 
and the calculated density matrix <r(t). Before estab- 
lishing such connection, we note that from a didactic 
point of view, it is easier to connect both methods 
described in this paper if we calculate the total field 
in the far-field zone, i.e., r ^> ri2, (ujba/c)^ 1 . We 
stress, however, that our conclusions are also valid 
for observation points r in the near-field zone. The 
total electric field operator is given in the far-field by [ljj 




(a) 




(b) 



FIG. 2: The coherent coupling (Figure a) and the incoher- 
ent coupling (Figure b) (in units of F) as a function of the 
dimensionless parameter x = uJbariz/c. We have chosen fi 
perpendicular to ri2, which is the case in the system studied 
in this paper. 



E(r,t) = E+(r,t) +E"(r,t) 



zM-) 2 E r^( d - 

47T£ V c / *-' n r- r. 



R 



i=l,2 

Sr(t-|r- ri |/c) 



+ H.c, 



(16) 



where the shorthand notation R, = r — has been used. 
Taking the expectation value and Fourier transforming 
yields 



III. THE T-MATRIX APPROACH TO 
SCATTERING OF LIGHT 



= {—) E^r^( d ^-^— ^— )x 

4=1,2 111 » 
e if |r-r<| /•+«> , _ 

-^-L *«" rt ((fir(*)>+H.c.). 

(17) 

If we now write (using the Green's funtion G defined 

as 

(d afc - R, (R ' ff2 dab) ) - -47r|r|d ai ,e-^l r l < G (r - r,, w) • /z, 

**4 

(18) 

then 



We now reconsider the system shown in Figure ^ We 
model both atoms as point-dipoles with dipole moments 
given by (JTJ, and study the dynamics of the system 
using a T-matrix approach instead of a density matrix 
approach. The total field E at any position in the 
three-dimensional coordinate space can be written as 
the sum of the incident field and the field scattered by 
both atoms |20| : 

E[u>, r] = E q [uj, r] 

+ J J Go(w,r — s)T(oj, s,u)E [uj,u]dsdu, 

(22) 

where the integrations are over the full three-dimensional 
coordinate space. The 3x3 dyadic Green's function 

Go(ui,r) has as representation in coordinate space pi) : 



[E[rM) = -^(v) 2 £ (tfo(r-r i>W ). M )x 

i— 1,2 

"((5,3))+H.c). (19) 



e \ c 

+00 



In the steady-state regime, the coherences in the system 
(the off-diagonal elements of the density matrix) oscillate 
in phase with the incident field. We can easily study 
the behavior of the coherences by transforming to the 
interaction picture 



S+ = e- iuot Sf, 

sr = e +iuot sr, 



e{l,2}, 



(20) 



such that the expectation values (^S i j are time- 
independent in the steady-state regime (see, e.g., 
The field component at lj — loq (the frequency of the 
incident field) is then given by 



E[r,w ] = (E[r,w 



4=1,2 



(21) 



Expression (|21|l expresses the total field E on resonance 
as a function of the scattering properties of the atoms, 
without making any assumptions regarding the indepen- 
dence of both scatterers. 



G ((j,r) 



.(/(3) 



uj 2 /c 



1 V ® v) e 



4irr 



iLfpu^rUV) +P(i-r)r2>r/t*), (23) 
rr V c c / 



with 



1 1 



P(z) = l-- + — and P(z) = -l + ? , 

z z* z z z 

(24) 

where ® denotes the tensor product. The 3x3 unit 
tensor is denoted as 1^ . We note that the Green's 
(|23|l function is related through l|24|) to the coherent 
coupling <!)( 12 ) and incoherent coupling r*- 12 ) which arise 
in a master equation approach (see expressions (|1 5|) ^ by 



/x • G (uj,r 1 - r 2 ) • /i 



1 u ba /c 



6?r r 



bo 



^(12)_ ir (12)\ 



(25) 



indicating that coupling through the vacuum is inti- 
mately connected to ordinary free-space propagation. 
The T-matrix T of the total scattering system is given by 

T{u,r,r') = {r\T(u)\r') 

= H®fit(ui) 5(r — ri)S(r' — rj)x 

([/W-B^n.raJtH]- 1 ^-, (26) 



where the delta-functions clearly depict the local char- 
acter of both scatterers. The theJ2 x 2 unit tensor is 
denoted by 1^ . The 2x2 matrix D is given by 



D(uj,r 1 ,r 2 ) = M • G (uJ,r 1 ,r 2 ) ■ fi 



1 

1 



(27) 



Each atom has a T-matrix element t(uj) describing its 
scattering properties. In the system we are considering 
the T-matrix elements of both atoms are equal because 
the atoms are identical. In the limit of low-intensity 
incident fields, the T-matrix element t(u>) is given by the 
well-known expression for the dynamic polarizability of 
a two-level atom 221 



3tt 



UJ - UJ ba 



(28) 



For incident waves of higher intensities, both atoms will 
show saturation effects. The extension of expression 
(12*51 beyond the low-intensity limit is given by [23|, |24( 



3tt 



i 



UJba I C 



(29) 



which takes into account the loss of photons at frequency 
uj due to inelastic scattering; the inelastically scattered 
photons at frequencies uj' ^ uj are not described. The 
saturation parameter s is defined as (see, e.g., 0) 



n 2 /2 



{uj ba -ujy + T*/^ 



(30) 



with the Rabi frequency = \d ab ■ Eq\/H. With the 
T-matrix element l|28l) we can rewrite the total field Ij22(l 
as 



E[oj, r] =E [uj, r] 



t(uj) 



1 - t{w)n ■ Gq(uj, n,r 2 ) ■ A* 



^2 G a (uj,r - n) ■ n® A* ■ E [u},rj]. (31) 

Equation (|31|l is the result of a T-matrix calculation, ex- 
pressing the total field E as a function of the incident field 
Eq and the scattering properties of both atoms, explicitly 
assuming that they are independent from one another. 



IV. COMPARISON OF BOTH METHODS AND 
DISCUSSION 

We now compare expressions <|21(1 and l|31|) and discuss 
their discrepancies. We consider three regimes. 



A. Low-saturation incident field 

If the saturation is very small (s <C 1), the steady- 
state value of the density matrix can be calculated 
analytically. The resulting steady-state expectation 

value (S^, i G {1,2} is, using the linear expression 



lim i 



s; 



t(uj) 



67rr ba i _ t(u>)n ■ G q (uj, n, r 2 ) ■ H 



(32) 



which yields, if substituted in l|21|) . the same expression 
()31J) as given by a T-matrix approach, hereby proving the 
validity of the T-matrix formalism for low-intensity inci- 
dent fields for arbitrary separation between both atoms. 



B. Atoms in each other's far-field: rrz 3> (uJba/c) 1 

In the limit of large separation between both atoms 
r 12 3> (uJba/c)^ 1 , both the coherent coupling 5^ 12 ^ and 
the incoherent coupling r^ 12 ' reduce to zero. In this 
regime, the density matrix a obeys 



Tri(o-) = Tr 2 ((r) = cr , 



(33) 



where <7o is the (2 x 2) density matrix corresponding to a 
single two- level system (see, e.g., and Tri stands for 
the trace over atom i. In other words, if both atoms are 
in each other's far-field, the single-atom density matrix 
factorizes out in two single-atom density matrices. This 
result corresponds with the intuitive expectation that the 
independent scattering approximation l|31|l is valid for 
large separation. 



C. Arbitrary separation and incident field in the 
saturated regime 

If the saturation parameter s is of order 1 or larger and 
the separation between the atoms is of the order of (or 
smaller than) (wb a /c) _1 , no easy analytic expressions for 
the density matrix exist. However, the density matrix 
and the total field can still be evaluated numerically. 
In Figure - El we have calculated fi ■ E[ujQ,r] using 
expressions (|21|) and l|31|) . Obviously, the numerical 
value of the total field component fi ■ E [ujq , r] depends on 
numerous parameters such as the incident field strength, 
the distance between both atoms, the evaluation point 
r and the frequency of the incident field. It is not 
our intention, however, to fully explore the entire 
(ujba — uJo,r,ri2,fl) parameter space. We only want to 
point out that differences between both methods do 
arise, and explain why. 



C Arbitrary separation and incident Held in the saturated regime 
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FIG. 3: The real part of the field component /x-25[o;o]/-E (0) on 
resonance, evaluated at the perpendicular bisector of ri2 • The 
Rabi frequency of the incident field is chosen fi/T = 1. The 
field has been calculated using both the T-matrix formalism 
(solid line) and the density matrix formalism (dotted line). 
The field has been plotted as a function of the detuning S, for 
(a) atoms in each other's far-field, and for (b) atoms in each 
other's near-field. 



FIG. 4: The imaginary part of the field component fj, ■ 
on resonance, evaluated at the perpendicular bi- 
sector of ri2. The Rabi frequency of the incident field is 
chosen Q/T = 1. The field has been calculated using both 
the T-matrix formalism (solid line) and the density matrix 
formalism (dotted line) . The field has been plotted as a func- 
tion of the detuning 5, for (a) atoms in each other's far-field, 
and for (b) atoms in each other's near-field. 



As an example, the field has been evaluated at the 
perpendicular bisector of r%2- The field is shown in units 
of 

E (o) _ Jtt . /g o(r _ ri) + g o(r _ A . ^ (34) 

UJ ba /C \ J 

Figure and 01 show the field as a function of the detun- 
ing 6 = 2(u!ba — wo). Figure [S] and show the field as a 
function of the Rabi frequency f2, expressing the inten- 
sity of the incident field. It is clear that if both atoms are 
in each other's far-field (r i2 ^> (wb a /c) _1 ), or if the satu- 
ration parameter is very small (s <C 1), the predictions of 
the T-matrix approach and the density matrix approach 
coincide, as explained above. However, we see that for 
atoms in each other's near-field (r 12 ~ (w(, a /c) _1 ), the 
predictions of both approaches deviate for increasing sat- 



uration parameter of the incident field. At very high sat- 
uration (s 1), both approaches again yield the same 
results. 

We can explain the difference between both approaches 
as follows. For low-intensity incident fields, both atoms 
scatter light purely elastically. At higher intensities, 
more light is scattered inelastically (and less light elas- 
tically). In a T-matrix approach, the coupling of both 
atoms with inelastically scattered light is not taken into 
account. This explains why, in the case of a single scat- 
terer, the T-matrix approach is exact, why for multiple 
scatterers, problems can arise when using T-matrices. 
The reason for this is that in the regime where the T- 
matrix approach breaks down, it is impossible to quan- 
tify the incident light on each scatterer, since the inci- 
dent light to a certain extend consists of incoherent light, 
generated by multiple inelastic scattering between both 



C Arbitrary separation and incident Held in the saturated regime 
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FIG. 5: The real part of the field component fi ■ E[uio]/E^ 
on resonance, evaluated at the perpendicular bisector of ri2. 
The detuning of the incident field is chosen 5 = 0. The field 
has been calculated using both the T-matrix formalism (solid 
line) and the density matrix formalism (dotted line). The 
field has been plotted as a function of the Rabi frequency Q, 
for (a) atoms in each other's far-field, and for (b) atoms in 
each other's near-field. 




(b) ri 2 =0.8K/c)- 1 

FIG. 6: The imaginary part of the field component fi ■ 
on resonance, evaluated at the perpendicular bi- 
sector of ri2. The detuning of the incident field is chosen 
5 = 0. The field has been calculated using both the T-matrix 
formalism (solid line) and the density matrix formalism (dot- 
ted line). The field has been plotted as a function of the Rabi 
frequency fi, for (a) atoms in each other's far-field, and for 
(b) atoms in each other's near-field. 



atoms. It is then no longer possible to define a single- 
scatterer scattering property, since coherences arise be- 
tween different scatterers, which are induced by inelas- 
tically scattered light. Finally, at very high intensities 
of the incident field, both atoms will be completely sat- 
urated, regardless of the approach one uses to calculate 
the field. This saturation explains why for large s, the 
results produced by both methods again coincide. 



V. EXTENSION 

We stress that the breakdown of the independent 
scattering approximation we describe here applies to 
many other systems. In general, all systems where the 
contribution of incoherent light is not negligible will 



experience a breakdown of the independent scattering 
approximation. If, for example, one wishes to describe 
atomic systems with gain l23l or the effect of saturation 
on coherent backscattering [23) HE- H3 using independent 
atoms, one is limited to large separations between the 
atoms. If one wishes to describe these systems for atoms 
in each other's near-field, then the atoms cannot be 
described (from a light scattering point of view) as 
individual entities; they are to be considered as one 
large indivisible system. Only a description taking into 
account all self-induced internal coherences, such as a 
density-matrix description, then produces correct results. 



VI. SUMMARY 

In this paper, we studied the scattering properties of a 
system consisting of two atoms in free space. The scatter- 
ing of incident monochromatic light has been considered, 
using both a saturable single-frequency elastic T-matrix 
approach and a density matrix approach. The former 
method assumes both atoms to be independent, while the 
latter method does not make such initial assumption. We 
found that in the regime of low saturation parameters, or 
for large distance between both atoms, both approaches 
yield the same predictions. However, if the atoms are 
in each other's near-field and the saturation parameter 



is high enough, the two approaches produce different re- 
sults. This discrepancy is due inelastic scattering of light, 
which is not taken into account in a T-matrix approach. 
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